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We compute the one loop self energy, in a locally de Sitter background, for a 
massless fermion which is Yukawa-coupled to a massless, minimally coupled scalar. 
We then solve the modified Dirac equation resulting from inclusion of the self energy. 
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, We find faster-than-exponential growth in the fermion wave function, consistent with 

' the production of fermions through a process in which a scalar and a fermion-anti- 

■ fermion pair are ripped out of the vacuum by inflation. 
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I. INTRODUCTION 

Consider quantum field theory in the presence of a homogeneous, isotropic and spatially 
flat geometry, 

ds^ = —dt^ + a^{t)dx ■ dx . (1) 

By virtue of spatial translation invariance we can label particle states by their constant wave 
vectors, k. However, the physical 3-momentum associated with this label is k/a{t). And if 
such a particle has mass m, its energy is. 



E{t,k) = \/m^ + \\k\\ya^{t) . (2) 



* At the end of the paper an addendum is appended which contains Erratum for this paper. 
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Now consider the appearance, at time t, of a pair of virtual particles with wave vectors 
±k. This process conserves 3-momentum but not energy. According to the energy-time 
uncertainty principle this violation of energy conservation will not be detectable provided 
the pair disappears back into the vacuum before a time interval At such that, 

dt'2E{t',k) <1 . (3) 

For any m 7^ the integral grows without bound, so At is finite for massive particles. 
However, if m = and the scale factor grows rapidly enough, the integral remains bounded 
even for At — > 00. This means that massless virtual particles need never recombine. 

"Rapidly enough" turns out to mean d{t) > 0, which is a type of expansion known as 
"inflation" . Whether or not inflation leads to a significant amount of particle production 
for a given species depends upon the rate at which virtual particles of that species emerge 
from the vacuum. It turns out that the emergence rate for particles which are conformally 
invariant falls like l/a{t) l^j. This means that any conformally invariant and massless virtual 
particles which happen to emerge from the vacuum will become real, but not many emerge. 
Hence there is negligible direct production of conformally invariant particles during inflation. 

Most familiar massless particles possess conformal invariance. In particular, photons 
are conformally invariant in D = 4 spacetime dimensions, and massless Dirac fermions are 
conformally invariant in any dimension. On the other hand, neither massless, minimally 
coupled scalars nor gravitons are conformally invariant, and both are produced copiously 
during inflation. This effect is observable; the production of nearly massless, minimally 
coupled scalars is believed to be responsible for the primordial anisotropies in the cosmic 
ray microwave background 0, 0] . 

One naturally wonders what happens when a massless but conformally invariant particle 
interacts with a massless particle which is not conformally invariant. The first step in 
answering this question is to compute the one particle irreducible (IPI) 2-point function 
of the conformally invariant field. One next uses this IPI 2-point function to correct the 
linearized equations of motion. The effect of inflation on the conformally invariant particle 
can be inferred by solving these equations for the plane wave mode functions. 

We have previously carried out such a study for photons interacting with a charged, 
massless and minimally coupled scalar 0, 0]. The result is that the mode functions of 
super-horizon wave vectors oscillate around a nonzero constant with increasing frequency 
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and decreasing amplitude This implies zero photon production and an enormous en- 
hancement in the 0-point energy, roughly similar to what happens for massive photons in 
Proca electrodynamics. The physical interpretation is that the inflationary production of 
charged scalars results in a diffuse plasma which inhibits the propagation of light. 

The purpose of this paper is to study what happens when Dirac fermions are made to 
Yukawa-interact with a real, massless and minimally coupled scalar. We will show that 
although the fermion mass remains zero, the mode functions for super-horizon wave vectors 
grow in a manner that is consistent with the production of fermions. The physical interpre- 
tation seems to be that inflation alters the constraint of energy conservation to permit the 
spontaneous appearance of a scalar and a fermion-anti-fermion pair. 

In Section II we compute the one loop fermion self energy. In Section III this is used 
to correct the Dirac equation. Although this equation is nonlocal, an asymptotic form is 
derived for the late time behavior of a spatial plane wave. In Section IV we show that the 
particle production inferred from the mode functions is consistent with a 3-particle creation 
process. Our results are summarized and discussed in Section V. 

II. THE ONE LOOP SELF ENERGY 

We begin by reviewing the conventions appropriate to Dirac fields in a nontrivial geom- 
etry. In order to facilitate dimensional regularization we make no assumption about the 
spacetime dimension D. The gamma matrices (6 = 0, 1, . . . , D — 1) anti-commute in the 
usual way, {7*, 7^^} = —2i]^'^I. One interpolates between local Lorentz indices [b, c,d, . . .) 
and vector indices (lower case Greek letters) with the vierbein field, e^b(x). The metric is 
obtained by contracting two vierbeins with the Minkowski metric, g^u{x) = efj_b{x)eyc{,x)rf"^ . 
The vierbein's vector index is raised and lowered by the metric (e^^ = g'^'^e^b) while the local 
Lorentz index is raised and lowered with the Minkowski metric (e^^ = rf"^e^,^. The spin 
connection and the Lorentz representation matrices are, 

A,bc ^ e\ (e.,,, - TP^^e,:) , J'^ ^ [7^ 7I • (4) 

Let represent a real scalar field and let ipi{x) stand for a Dirac field. In a general 
background metric the Lagrangian we wish to study would be. 
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where i/j = ip'^j^ is the usual Dirac adjoint and / is the Yukawa couphng constant. The 
geometry of interest is the very special form associated with the homogeneous and isotropic 
element (^. By defining a new time coordinate drj = dt/a(t), the metric of this geometry 
can be made conformal to the Minkowski metric, 

ds"^ = a^(^—drf + dx ■ dxj . (6) 

A convenient choice for the associated vierbein is, e^t = ar^^b. With these simplifications the 
spin connection assumes the form, 

e^b = aT]i,b =^ Af^cd = [jli^cdd - r]f,ddc) ln(a) . (7) 

And our Lagrangian reduces to, 

Cconi = ~a'^^d^(t>du(t>r]^'' + (a'^ii)i^^d^{a'^i?) - fa'^H^i; . (8) 

Having defined the Lagrangian we now give the propagators and vertices which comprise 
the position space Feynman rules. Of course vertices are always straightforward. The (pipiipj 
vertex and the fermion field strength renormalization are, 

— ifa^Sij and i6Z2{aa)^~'yijid^6^{x — x'), (9) 

where a = a{t) and a' = a{t') are the scale factors evaluated aX x^ = f] and = 77' , 
respectively. We do not need scalar mass and field strength renormalization. 

The fermion propagator is best expressed in terms of the conformal coordinate interval, 

Ax'^{x; x') = \\x — x\\'^ — {\i]—r]'\ — i5Y . (10) 

Note that we label the position in spacetime through the D-vector x^ = {1], x). We see from 
the second term of (jHI) that, for / = 0, the combination a^~ilj behaves like a free, massless 
Dirac field in fiat space. It follows that the propagator of ip is just a conformal rescaling of 
the fiat space result, 

[ 4:71- L J J 2n- [Ax\x;x')]- 

(11) 

Here Ax^ = rj^^^x^ — x'^). The split index notation in i[iSj\{x\x') indicates that the first 
index {%) transforms according to the local Lorentz group at the first coordinate argument 
{x^) whereas the second index (j) transforms at the second argument {x'^). 



iSj 
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It has not been necessary so far to make any assumption about the form of the scale 
factor. However, the free scalar is not conformally invariant, which makes its plane wave 
mode functions depend upon a{t) in a complicated manner. Although the solution for 
general a(t) has recently been obtained the Fourier mode sum has only been worked out 
for certain cases. We therefore specialize to locally de Sitter inflation, 

1 



ait) 



Hr] 



(12) 



In this geometry there is a simple relation between i{x; x'), the invariant length from to 
x'^, and the conformal coordinate interval Ax^(x;x'), 



4 sin^ -Hi{x] x') = aa'H^Ax^{x; x') = y{x; x') 



(13) 

We refer to y{x; x') as the de Sitter length function. The scalar propagator can be expressed 
in terms oiy = y{x; x') and the two scale factors, 



iA(x: x') 
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(14) 



n=l V 2 

Note that the homogeneous terms (i.e., y") and the constant terms have slightly different 



proportionality constants from our previous expression |^ . This has been done to make 
()14|) valid for regulating a spacetime in which D will ultimately be taken to some dimension 
other than four y^j. In D = 4, the propagator reduces to a more familiar, elementary 
form. 



iA(x; X 



\n{H^Ax^' 



1 
4 



ln(2) 



(15) 



47r2 I Ax2 2 

The diagrammatic representation of the one loop fermion self energy is given in Fig. 
Computing it is simple: just multiply the various vertices and propagators. 



{x]x) = {^ifa^ 5ik^i kSe_ {x;x')(—ifa^6£j^iA{x;x')+i6Z2{aa) 2 ■y-'jid^6^{x—x') . 

(16) 

To renormalize in four dimensions, set D = 4 — e and then expand in e. Only terms which 
are not integrable over d^x' need to have e arbitrary. 



i\D- 



iSj 



fx; x')iA(x: x') 
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2%4 4 



2^ {Axy 

U^H^Ax^ 
V 4 



i-fijAx^ 
(Ax')2 



+ 0(e). (17) 
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FIG. 1: T/ie one /oop fermion self energy. Fcrmion lines have an arrow, scalar lines do not. The 
final diagram gives the contribution of field strength renormalization. 



The next step is writing the factor of Axn as a derivative, 



)■ 



-Pi 



1- /Kf^^Ax^l^ 



(18) 
(19) 



One then expresses inverse powers of Ax^ as powers of 9^ = rj^'^d^di, acting on logarithms, 

Ax2 8 W ^ ^ 

It is at this stage that the renormalization scale /i appears and that the ultraviolet divergence 
is segregated to a local term, 

1 1 



Aa;2 



---9^ In2(/x2^a;2)-21n(/x2^a;2) 



efl -e)r(l - 



6^{x-x') + 0{e) 



(21) 



Defining ■JijO/^ =^ij and combining terms we have. 



-f{aa' 



iSj 



ix;x')tAix;x')^i^^^t^.,d' 



ln"(/x"A,T") -21n(/iW) 



pH\aa')i .. 



^..a2ln2(Yi^'Aa;2) + 



?/V^(aa0^r(l-f) 
287,4 ^rzj- - V 4 y ■ e(l-e) 

The simplest renormalization condition is, 

/V^ r(i - f ) 



i^ijS'^ix - x') + 0(e) . (22) 



SZo 



2%2 



6(1 -e) ■ 



(23) 



We obtain the renormalized self energy by taking the limit e — > 0, 

Piaa'y 



v->rcn 



{x; x') 



10 4 ^^J^Un\^^'Ax')-2H^/Ax') 



P{aa') 



In(aa') i^ij6^{x — x') — 



pH\aa' 



■^d'ln^^H'Ax')+0{p) . (24) 



2%2 -^--^^y^J-^- 287r4 r.j V4 

The first term is the conformally rescaled flat space result. The second term is the well- 
known contribution from the conformal anomaly. The intrinsic de Sitter result is the third 
term. It is distinguished from the other two by its extra factor of aa' . 
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III. SOLVING THE MODIFIED DIRAC EQUATION 



The self energy modifies the Dirac equation to include quantum corrections. The modified 
equation takes the form, 



3 

an 



X 



(x; x')'iljj{x') = . 



(25) 



It describes how the mode functions change due to quantum corrections. From these mode 
functions one can infer particle production, and also kinematical properties such as mass. 

There is some ambiguity about what it means to "include quantum corrections," and 
each of the possible meanings corresponds to employing a different self energy functional in 
the modified Dirac equation ()25|) . The in-out self energy corresponds to a particle which 
is surrounded by free vacuum in the distant past, and for which quantum effects conspire 
to leave it surrounded by free vacuum in the distant future. This self energy is appropriate 
to fiat space scattering experiments, but it is not very realistic for cosmological settings. A 
better choice for cosmology is to release the particle surrounded by free vacuum at some 
finite time and then let it evolve as it will. For that situation it is the Schwinger-Keldysh 
self energy which is appropriate [l^ . 

At one loop order the Schwinger-Keldysh self energy is simple to construct from the in- 



out self energy ()24|) . We will give the technique without deriving it 



Each vertex point 



is endowed with a polarity ±. In place of the conformal interval (fTIH) we define. 



Axl+ = Ar^ - {\Ari\-i5y 



(26) 



where Ar = ||x — and At] = r]—r]'. At one loop order we do not require the — h or 
variations. The ++ and H — self energies are. 



(x+; x'_^) ■ 
P{aa' 



2i07r4 



^^\n{aa')i^,,5\x-x') - ^ ^J^^ ^ 



^dHnH^H'Axl 



+ 0(r), (27) 



P{aa' 



2i07r4 



^ijd^ \ii\n^Axl__) - 21n(/i2Ax^ 



fH\aa')l 
2%4 



^.,dHn'{y^H'Axl_)+0{f). (28) 



The Schwinger-Keldysh self energy is the sum of these. 



SK 



(x+;x+) + 



[x+;x_ 



(29) 
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Our explicit result for [jS^ is 



hjd^[e{A7])e{A7]-Ar)\n[e^y{A7]'^-Ar'^)]}- 
pH\aa') 



Pjaa') 
257r2 



In (aa') i^ij5^{x — x') 



i^ijd^i^9{Ari)e{Ari-Ar) \n[^H^{A7f - Ar 



0{f). (30) 



Note that it vanishes whenever the point x'^ is outside the past light-cone of x^. This means 
that although the modified Dirac equation ()25|1 is nonlocal, it is nonetheless causal. Another 
important property of [jS^^](x;x') is that each term contains only a single gamma matrix. 
This means that it conserves helicity, and hence that the fermion remains massless. 

Note also that the first term of ()3U|). the flat space part, has the same number of scale 
factors as the tree order differential operator in (j25|) . Perturbation theory only makes sense 
if the coupling constant is small, so we may assume p <^1. This means that the flat space 
part of the one loop self energy can be ignored. So too can the conformal anomaly, which 
is only enhanced by a factor of ln(a). However, we cannot necessarily ignore the de Sitter 
correction. Although it is also down by p, it is enhanced by a scale factor, which becomes 
enormously large during inflation. It turns out that higher loop corrections — which are 
down by more powers of / — can never give more than this one factor of a. It therefore 
makes sense to keep only the one loop de Sitter contribution of (fHUjl . 



pH^{aa')l 
¥7^ 



i^ijd'^{e{Ar])e{Ar]-Ar) In '^H'^{A7f-AP 



(31) 



This is the "self energy" we will henceforth employ in the modified Dirac equation (j2| 

Because the background is spatially translation invariant it makes sense to look for plane 
wave solutions, 

^i{r],x) = a-^^Xiijl.k)e'^'^ , 
To take advantage of helicity conservation we employ chiral gamma matrices, 

h::) • -(:;■)■ 

and we factor out the tree rj dependence of the 2-component helicity eigenstates, 

/X^(r/,A;)e±*'=n 



(32) 



(33) 



with 



For the left-handed spinor the result of performing the spatial integrations is. 



PH 



doXLiv,k) + I dr]'axL{v',k)e 



\n{HA7]) + 



3 1 
4^2 Jo 



2k At] Q±iT_'^ 

dr 



(34) 



0. (35) 
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Here r]i = —1/H is the initial conformal time, which corresponds to t = 0. The analogous 
equation for the right-handed spinor with ± helicity is, 



2ln{HAr]) + -- 



dr 



0. (36) 



Recovering the full time evolution of nonlocal equations such as ()35II36|1 requires numerical 
integration. However, it is straightforward to get the asymptotic behavior for late times. 
First change the independent variable to the number of e-foldings, = ln(a). Now redefine 
the dependent variable as hf^ ^{N,w), where w = k/H. The equations which result are. 



f 2 

dNh^{N,w) + (^) / dN'h^{N',w)exp\Ti2w{e-^'-e-^) 



'27r/ Jo 
X < In 



3 1 r2wie-'^'-e-^) 



+ 4 + 2 



dr 



. 



(37) 



The equations for the right-handed spinor are obtained by replacing =pz2w with ±2iw and 
±ir with =l=ir. 

To obtain the asymptotic form for large we simplify fl37|) by dropping terms in the 
integrand which are exponentially small in for fixed A^' < A^, 

d^ht{N,w) + (^) dN'ht{N',w)t 

^ /TT / .In 



^2n^ Jo 

At large A^ we see that the system has the form 



A 2 Jo 



dr 



T 



. (38) 



N 



dNhUN,w)= / dN'gUN',w)h^{N',w) 



(39) 



where the function goo{N,w) is. 



3 1 /-a^-e-^ 

dr 



4 2 Jo 



(40) 

Because the function goo{N,w) grows, the integral is dominated by its upper limit and we 
can obtain the asymptotic form of hoo{N,w) (up to a proportionality constant) using the 
WKB technique, 

(AT, w) g-J (AT, w) exp [ fdN'^g^N', w)] ^ (a^ - ^)"^ exp[£(Ar - ^)^]. (41) 

Each of the various helicities and polarizations, hf^{N,w) and h^{N,w)^ approaches a con- 
stant times this same limiting form hoo{N,w). 
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IV. THE PHYSICS BEHIND OUR RESULT 

A. Massless minimally coupled scalars 

Most familiar particles become (classically) conformally invariant when their masses are 
taken to zero. The exceptions are gravitons and massless, minimally coupled scalars. The 
Lagrangian density for the later is, 

>Cmmcs = -ld^<t>d.(t>9'"'V^ = —a^d^cPd^cPr)^" . (42) 

The Lagrangian is the spatial integral of the Lagrangian density. Using Parseval's theorem 
to convert to Fourier space, and expressing the field in terms of co-moving (rather than 
conformal) time, gives a form we can recognize, 

i^MMCS^ / ci^a;£MMCs = ^/^{a='(t)|^(t,plP-a(t)lbiri0(t,p)P}. (43) 

Each wave vector p represents an independent harmonic oscillator with a time dependent 
mass, ■m{t) ~ a^{t), and frequency, u!{t) — \\p\\/a{t). 

For the case of de Sitter inflation (a{t) — e^*) the time dependence of the part of the 

— * 

field with wave vector k takes the simple form, 
Ht,P) = f{t,p)ai^ + nt.pW{-p) where f{t,p) = " ^Y'""'- (^4) 

To understand the meaning of the canonically normalized creation and annihilation operators 
operators a{p) and Q!^(p), note that the minimum energy in wave vector p at time t is 
\{h)uj{t). However, the state with this energy does not evolve onto itself. Indeed, this 
theory has no stationary states! A reasonable "vacuum" is the state that was minimum 
energy in the distant past. This is known as Bunch-Davies vacuum and it is defined by 
a(p)|n) = 0. 

For Bunch-Davies vacuum the 0-point energy in wave vector p is, 

E,{t,p) = \a\tmt,p)\' + \amp\?\mp)? = + ■ (45) 

The first term on the far right represents the irreducible, minimum energy. The second term 
gives the extra energy due to inflationary particle production. Since the energy of a particle 
of wave number p is ||p||/a(t), we can easily compute the number of particles, 

- mj . (46) 



11 



As one might expect, this is much less than one at very early times, and it becomes of order 
one when the wave number just begins to satisfy the uncertainty bound 

B. Production of fermions 

Although we have seen that inflation produces enormous numbers of massless, minimally 
coupled scalars, the conformal invariance of free Dirac theory implies that there can be no 
comparable, direct production of fermions. However, it is still possible to make fermions 
during inflation by allowing them to interact with a massless, minimally coupled scalar. To 
see this we begin by giving the free field expansions of the interaction picture fields, 

Mt,x) = J-0^,{a{p)f{t,p)e'P^-'^ + a\p)r{t,p)e-'^^-'^^ (47) 
Mt, S)= / E r)^e-/^'^+--^ + y (g, r)^e--/^"---j , (48) 

V^,(t, f) = / ^ E |7(^^, sf-^e^^'^'^^^^-' + PKk, .)^e-^/^-'^^-j . (49) 

Note that we have absorbed a factor of at into ipj and ipj so that they agree exactly with the 
flat space expansions in conformal coordinates. The spinor wave functions obey the usual 
spin sum formulae of massless fermions in flat space, 

E u{k, s)u{k, s) =/if = E "^(^5 s)v{k, s) . (50) 

s s 

The nonzero relations of the operator algebra are, 

'a{p),a\p')] = {27i)H\p-p'), (51) 
{P{k,s),p\k\s')] = {2iif5\k-k')5,, = {^ik,s),^^ik',s')} . (52) 

We are studying the same Yukawa theory as in sections II and III and, as before, we 
release the universe in free Bunch-Davies vacuum at t = 0, 

aip)\n) = = P{k, s)\n) = 7(g,r)|^]) . (53) 

Therefore the operator which maps the various free fields into interacting fields is, 

U{t) = T jexp [-ifj^dt'jd^x(f)i{t', x)'^j{t', x)%l)i{t' , x)] | . (54) 
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Although Heisenberg states such as \Q) do not change with time, what these states mean is 
determined by operators, which do change with time. In particular, the mapping. 



J{t)^{t,x) = U\t)^i(t,x)U{t) 



suggests a reasonable definition for the operator which annihilates fermions at time t. 



(55) 



f3{t,k,s) = U^{t)P{k,s)U{t) . 



(56) 



We therefore arrive at an expression for the probability that \Q) contains a fermion of wave 
number k and spin s at time t, 



Ppit, k, s) {2Tif5\k - q)6rs = (3\t, g, r)(3{t, k, s) 



Q) = (Q 



U^it)P^iq,r)P{k,s)U{t)\n) . 

(57) 

Expression (|H7j) is straightforward to evaluate to lowest order. One first commutes the 
creation and annihilation operators through the factors of f/^(t) and U{t), 



Ppit,k, s)i27r fS%k-q)S, 



mr).U{t) Pik,s),Uit) 



n) . 



(58) 



From the definition of U {t) , and the free field expansions, we find. 



)3{k,s),Uit)]=^ 

^ V2k Jo 

The free field expansions imply. 



0/(^", y)i'i{t'\ y)u{q, r)(f)i{t', x)u{k, s)i)i{t', x) 

dH 1 



d^p 
(2^ 



nt',p)f{t",p)e 



X 



,0 ^3 Of E ^(^"i rMk, S)V{1 s')e-^^/Ha'+U/Ha"^Ui^-y) ^(gg) 



The integrals over x and y enforce q = k and i = —{k + p). With these identifications we 
can set r = s and spin average to eliminate the fermion wave functions. 



i Y: v{l s')u{k, sMk, s)vi£, s') = ^TT[f }t 



Putting everything together gives. 



\k+p\\k + {k + p)-k 



(61) 



P,{t,k,s) = ^^fdt'fdt''(4^J\t\p)f{t\p)e-^^lh~ltrrM^^^^^^^^ 



2k 



JO 



(27r)= 



k^ 



[k + p) ■ k 



\k+p\\ 



+ 0(f). 
(62) 
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The time dependence of the integrand of ()62p 



tp 



tp 



1 

' Ha" 



||p+A:||+p+fc) 



(63) 



indicates that the lowest order effect is driven by the physical process in which a scalar and 
a fermion-anti-fermion pair emerge from the vacuum. In fiat space quantum field theory this 
process would be forbidden by energy conservation. It is crucial to note that this is not put 
in by hand: even in fiat space one gets an expression analogous to ()62|) . but it goes to zero 
as the time t tends to infinity. What drives the fiat space expression to zero is the fact that 
no value of the scalar wave vector p gives zero for the total fiat space energy, +p + k. 

Hence the temporal integrations result in rapid oscillations that tend to cancel. As t grows, 
the result approaches an energy delta function that cannot be saturated. In our locally de 
Sitter background the inverse scale factors redshift the physical energy to zero so rapidly 
that there are only a finite number of oscillations. 

Note the essential role played by the break;ing of conformal invariance. Since free fermions 
are conformally invariant, and the Yukawa interaction is not inconsistent with conformal 
invariance, the effect of conformal breaking is concentrated in the scalar wave function, 
fit,p). Had we instead used a conformally coupled scalar the integrand would go to. 



1 



-H75^-lj^)(l|p+fc||+P+fc) 



(64) 



2a'a"p 

In this case the inverse scale factors suppress the result. 

To give a good approximate computation of Ppit, k, s) we note, from the preceding dis- 
cussion, that scalar wave numbers with p > HA = Hmin{a', a") engender oscillations in the 
energy exponent. We can therefore cut the scalar wave number off above this limit. (The 
infrared cutoff of p > H derives from the finite coordinate range of our spatial manifold 
j^.) The same energy exponent also gives rise to oscillations when k > Ha' or k > Ha", 
which cuts off the lower limits of the temporal integrations at t',t" >T = ln{w)/H, where 
w = k/H. With these cutoffs in place we can neglect terms that fall off like 1/a' or 1/a", 



f 



2 rt 



Pj3(t,k,s) — > — dtj dt" I dpp — 7 



k Jt 



HA 



H^ 



H 



d^p 



2p3 J (27r 



The angular integration is straightforward. 



r d'^p 
(2^ 



k + 



{k + p) ■ k 

\\k + p\\ 



27r2 



9{k-p) 



p^ - 
3P- 



{k + p) ■ k 

\\k + p\\ 



+ e{p - k) 



+ 0{f 



(65) 



1 + 



2ki 
3p- 



(66) 
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Now change variables from time to e-foldings, and take advantage of the symmetry between 
t' and t". 



Pp{t, k, s) 



f2 pN 

■' ' dN- 



27r27ln 



7 dN" ^ 

■/Into Jh P [ 



9{k-p) 



3k^ 



+9{p-k) 



^ 2k 
3p 



+ 0{f). (67) 



Performing the integrations and retaining only the leading terms for ^ \\i{w) ^ 1, we 
obtain, 



P (t, k,s) ^ dN' dN"{N" + Hw) + - + ...} + Oif) 



Utt- 



;{N' + ...}+0{f 



(68) 
(69) 



As one might have expected, our leading order result for P/3(t, k, s) is independent of s and 
w. Upon taking a square root of Ppit^k^s), one gets a quite good approximation for the 
exponent of the WKB mode function (j4ip we found in Section III. 



V. DISCUSSION 



We have found that the effect of inflation on Yukawa-coupled, massless fermions is in some 
ways opposite to the effect on the photons of massless scalar QED. Whereas the photon 
becomes massive, but suffers no significant particle production [7J, we find that fermions 
remain massless but experience enormous particle production. The physical mechanism is 
simple: inflation alters the constraint of energy conservation so that three massless particles 
can spontaneously emerge from the vacuum. The Yukawa interaction mediates this for a 
scalar and a fermion-anti-fermion pair, and the process has an appreciable amplitude when 
the scalar is minimally coupled. 

The mathematical analysis for Yukawa-coupled fermions is quite different from that re- 
quired for photons in scalar QED. Although the mode equations of both systems are nonlocal, 
the integrals of the fermion system are dominated by their upper limits because the fermion 
mode function grows. This means that we need only drop exponentially suppressed terms 
and then differentiate to extract a local equation for the asymptotic form. The photon mode 
equations were much more subtle because the photon mode function approaches a nonzero 
constant — as do classical photons — and the asymptotic effect of vacuum polarization 
shows up in the increasingly rapid oscillations around this value. 
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Of course fermions obey the Pauli exclusion principle, so at most one particle can be put 
in each momentum and spin state. The amplitude of the wave function cannot really grow 
past that point. Therefore the rate of production for any particular momentum and spin 
state must go to zero as the probability of this state being occupied approaches unity. This 
phenomenon is known as Pauli blocking. We do not see Pauli blocking because we are using 
only the one loop self energy. We expect Pauli blocking to appear at two loop order. 
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Addendum: Erratum for "Production of Massless Fermions during Inflation" 

Abstract 

We correct a sign error in the effective Dirac equation solved in Ref.^^. We discuss the 
changes this implies and their physical significance. 

In Ref. we calculated the one-loop self-energy for massless fermions which couple via 
a Yukawa interaction term to a massless, minimally coupled scalar during de Sitter inflation. 
However, within the system of conventions we used, the self-energy should contribute op- 
positely to the local derivative term in the quantum- corrected Dirac equation. The correct 
equation is, 

a^^^i^ia^^^ipix)) - J d^x'Jlretix, x')^(x') = . (70) 

Here a = —l/{Hri) denotes the scale factor during de Sitter inflation, Sjet is the retarded 
self-energy of the Schwinger-Keldysh formalism, and ip is the fermionic wave function. This 
is how equation (29) in |l2] should have read. 

The correct sign makes a simple and profound change in our analysis. The conformally 
rescaled, chiral components of spatial plane wave mode functions obey the equation. 



f 2 I 



■ 3 1 r2we-'^' e^"-l 

N' / dr 



A 2 Jo T 



0. (71) 
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where w = k/H, N = ln(a), / denotes the Yukawa couphng, and the approximation is 
that some exponentially small terms have been neglected. This should replace equation (38) 
in In the late time limit of large N = ln(a), both helicity components approach the 
same declining oscillatory form, 

>.u"-{N-?,r'\.,(±-MN-?,r). (72) 
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This replaces equation (41) in 

The correct asymptotic solution ()72j) stands in sharp contrast to the faster-than- 
exponential growth which pertains with the sign error! The physical interpretation of the 
corrected result seems to be that the inflationary production of massless, minimally coupled 
scalars drives the scalar field strength away from zero, which affects super-horizon fermions 
like a mass term. For a more complete study of the mass generation mechanism in the case 
when the scalar field is light and nearly minimally coupled we refer to Ref. jisl . 

Declining oscillatory behavior very similar to ()72j) is observed when studying late time 
dynamics of super-Hubble photons |3] . So the sign correction brings our results for Yukawa 
into pleasant conformity with what charged scalars do to photons during inflation 3,0,0,0]. 
In each case inflationary particle production forces the scalar away from zero, whereupon 
the effect on other fields is qualitatively similar (but not in every detail!) to what happens 
for nonzero scalar backgrounds in symmetry breaking. 

Note that the sign change reported here does not affect either the one loop Yukawa scalar 

n 

self- mass-squared or its impact on plane wave scalar mode functions [15]. So the conclusion 
of that study is still that the scalar cannot develop a large enough mass rapidly enough to 
affect what it does to the fermion. Nor does the sign change reported here affect the recent 
computation of quantum gravitational corrections to the one loop self-energy of massless 
fermions in Dirac + Einstein ||16|. Of course our sign change very much does affect the 
impact these corrections have on fermion mode functions [17]! 
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